Probing the absolute density of the Earth's core using a vertical neutrino beam 



Walter Winter 

School of Natural Sciences, Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540 

(Dated: February 2, 2008) 



We demonstrate that one could measure the absolute matter density of the Earth's core with a 
vertical neutrino factory baseline at the per cent level for sin 2 28is > 0.01, where we include all 
correlations with the oscillation parameters in the analysis. We discuss the geographical feasibility 
of such an approach, and illustrate how the results change as a function of the detector location. 
We point out the complementarity to geophysics. 
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Neutrino oscillation physics has entered the age of pre- 
cision physics, which means that the leading atmospheric 
and solar oscillation parameters are known to high pre- 
cisions and the next generation of long-baseline exper- 
iments will be highly sensitive to sub-leading effects. 
These long-baseline experiments send an artificially pro- 
duced neutrino beam of energy E on a straight baseline 
(length L) through the Earth to a detector, which is, de- 
pending on the neutrino energy, several hundred to many 
thousand kilometers away. In particular, the future po- 
tential high precision instrument "neutrino factory" Q 
leads to typical baselines ~ 700 km — 7 500 km relevant 
for neutrino oscillation physics. It is an interesting fea- 
ture of neutrino oscillations that the flavor conversion is 
sensitive to the electron density n e of Earth matter Q , 
which has been suggested to be used for neutrino oscil- 
lation tomography of the Earth's interior 0, 0, 01 • The 
electron density then translates into the matter density 
by p = n e m^/Y e with Y e = Z / A the "electron fraction" 
and mpf the nucleon mass. For a neutrino factory, the 
signal amplitude of such a measurement using the flavor 
conversion <-> v e is given by the parameter sin 2 2#i3, 
which has so far only been constrained by the CHOOZ 
experiment to sin 2 2#i3 < 0.1 [(J. Future long baseline 
experiments will find sin 2 2#i 3 > within the coming ten 
years if sin 2 2#i 3 > 0.01 (see, e.g.. 



Ref. 0). 



The most successful approach to the tomography of 
the Earth's interior has been seismic wave geophysics 
primarily using seismic waves from earthquakes to re- 
construct a profile of the Earth's interior. Most of the 
energy produced by an earthquake is deposited in shear 
waves (s- waves), which cannot penetrate into the Earth's 
(outer) liquid core (but might be partially converted into 
p- waves). A smaller fraction of energy goes into pressure 
waves (p- waves), which are propagated into the Earth's 
core, too. The waves are partially reflected at the mantle- 
outer core and outer core-inner core boundaries. There- 
fore, seismic waves are highly sensitive to density jumps 
and the positions of these boundaries. Seismic waves 
geophysics leads to a propagation velocity profile of the 
Earth's interior, which can be translated into a density 
profile with the equation of state. This conversion is 
based on a model for the shear/bulk modulus, which 




FIG. 1: (Color online) Positions of three of the major poten- 
tial neutrino factory laboratories, (typical) L = 3 000 km de- 
tector sites (dashed curves) , as well as potential detector sites 
with outer core crossing baselines (below thin solid curves), 
and inner core crossing baselines (within thick solid curves). 
The colors of the curves represent the different laboratories. 



implies uncertainties. In fact, the most direct informa- 
tion on the matter density distribution comes from the 
Earth's mass and its moment of inertia about the polar 
axis, which are, however, not uniquely determining it. 
In particular, there are many open questions about the 
Earth's inner core (see, e.g., Ref. ||). Therefore, a mea- 
surement of the absolute electron density of the Earth's 
core could provide very complementary information to 
geophysics. 

From the point of view of a neutrino factory, this would 
require a "vertical" baseline, i.e., a decay tunnel which is 
vertical. Since there are only a number major high energy 
laboratories which are candidates for a neutrino factory, 
all of these on the northern hemisphere, the geographical 
aspect is another important part of this problem. As it 
is illustrated in Fig.^ there are indeed possible detector 
locations (i.e., on land instead of in water) for many of 
the major laboratories with core crossing baselines. For 
example, a baseline from CERN to New Zealand crosses 
the inner core of the Earth. 

The geometry and the quantities of interest of the 
Earth tomography problem are illustrated in Fig. [2] 
Though neutrino oscillations arc, in principle, sensitive 
to structural features of the density profile, realistic ex- 
periments with only one baseline have only very limited 
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abilities for a detailed reconstruction because they are in- 
sensitive to structures much shorter than the oscillation 
length in matter £|. Therefore, it is important that the 
relevant independent parameters be identified. In partic- 
ular, treating the (average) outer and inner core densi- 
ties pc and pic as independent parameters (c/., Fig. 
implies strong correlations 0. However, these two quan- 
tities are not really independent, since the total mass of 
the Earth is extremely well known. In addition, we have 
very good knowledge on the Earth's mantle from seismic 
wave geophysics. If we furthermore assume that the po- 
sitions of the mantle-outer core and outer core-inner core 
boundaries are very well known from the reflections of 
seismic waves, we can reduce the number of parameters 
to one (either pc or pic)- Below, we will argue that it is 
always reasonable to measure the average density of the 
"innermost" shell a baselines crosses. 

Neutrino oscillations are, for the slowly enough varying 
matter density within each shell pi (mantle, outer core, 
or inner core), to a first approximation determined by the 
baseline averaged density pf = 1/L L p(l)dl. Here p(l) 
is the density along the baseline. On the other hand, 
the average density of the Earth is determined by the 
volume averaged density, which is, for each shell, given 
by pY — 4-7T /V fft 2 p(r)r 2 dr. Here r is the distance from 
the Earth's center and V — 4/3(i?2 — -ft?) 71 " the volume 
of the shell. This means that the density within each 
differential shell dr is weighted with r 2 . For the total 
mean density of the Earth we then have 
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The constraint Me = const, i.e., p v — const because of 
the (approximately) known volume, together with the as- 
sumption of a known/fixed pY f , leads by differentiation 
with respect to pY c to ApY = —(Ric/Rc) 3 &f>Yc — 
—0.04 Apic- Thus, a large change in the measurement 
in pYc can be compensated by a very small change in 
Pq because of the volume averaging. The same is, in 
principle, true for the outer core and mantle densities. 
Since the effect on neutrino oscillations is proportional 
to the baseline length within each shell and not to the 
volume, each shell contributes to the total averaged den- 
sity p L by a similar magnitude. It is therefore rea- 
sonable to treat the density of the innermost shell as 
baseline crosses as the parameter of interest and cor- 
rect the (better known) density of the next outer shell 
by Ap v . In the Reference Earth Model (REM) 0, 
p c ~ 11.2 g/cm 3 , p v c ~ 10.9 g/cm 3 , pf c ~ 12.9 g/cm 3 , 
and Pic — 12.8 g/cm 3 (L = 2Re), which means that 
any model dependence in the conversions between these 
quantities only enters as sub-leading effect in the mea- 
surement precisions. 

Based on the above discussion, use pf c = pic as pa- 
rameter for inner core crossing baselines, and p c = pc for 
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FIG. 2: (Color online) The geometry of the core tomography 
problem (not to scale). 



baselines which only cross the outer core (c/., Fig. EJ) ■ In 
fact, because of the actual sensitivity to the electron den- 
sity, we measure the effective density p cS — J2i Pi Y^/0.5 
averaged over the baseline, where the index i refers to the 
contained elements (for example, for iron Y e ~ 0.464). 
Since for most of the elements in the core Y e is close to 
0.5 and therefore very similar, errors in the composition 
translate into errors in p only as second order effect. Wc 
assume the matter densities to be constant within each 
shell, i.e., we ignore the matter profile effect, which cor- 
responds to measuring the absolute normalization of the 
Reference Earth Model (REM) profile. We use a com- 
plete simulation of the neutrino factory Nu Fact- 1 1 from 
Ref. [Hi simulated with the GLoBES software In 
particular, we include statistics, systematics, and (con- 
nected and disconnected) degenerate solutions with the 
oscillation parameters, i.e., we marginalize with respect 
to the oscillation parameters. This procedure is neces- 
sary to test if the matter effect sensitivity survives a 
realistic simulation taking into account the insufficient 
knowledge on the oscillation parameters. The neutrino 
factory uses muons with an energy of 50 GeV, 4MW tar- 
get power (5.3 • 10 20 useful muon decays per year), and 
a magnetized iron detector with a fiducial mass of 50 kt. 
We choose a symmetric operation with 4 yr in each po- 
larity. For the oscillation parameters, we use the cur- 
rent best-fit values Am^ = 2.5 • 10~ 3 eV 2 , sin 2 26> 23 = 1, 
Am^ = 8.2 • 10- 5 eV 2 , and sin 2 28 12 = 0.83 We 
only allow values for sin 2 26*i3 below the CHOOZ bound 
sin 26»i3 < 0.1 6] and choose Sqp = as well as a normal 
mass hierarchy as input parameters, where the results 
should hardly depend on the choices of Sep and the mass 
hierarchy. Furthermore, for the leading solar parameters, 
we impose external precisions of 10% on each Ato 2 ^ and 

In principle, it is enough to have one baseline to mea- 
sure the average density of the Earth's core. However, if 
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FIG. 3: (Color online) The measurement of pic as function 
of the true value of sin 2 26ig at the la, 2a, and 3a confidence 
levels (from light to dark shaded regions). For the baselines, 
L = 2 ■ Re combined with L — 3 000 km is used. The dashed 
curves correspond to fixing the oscillation parameters, i.e., to 
not taking into account correlations and degeneracies. 



a neutrino factory is built, it's main purpose will be the 
determination of Sep and other sub-leading effects, which 
means that (for a muon energy of 50 GeV) L = 3 000 km 
is a typical baseline for that purpose (c/., Fig. ^ for po- 
tential detector sites). Therefore, we show in Fig. the 
precision of the measurement of pic with L = 2 • Re 
as function of the critical parameter sin 2 2#i3 in combi- 
nation with L = 3 000 km to reduce the errors of the 
oscillation parameters. Note that we have chosen the 
longest possible core-crossing baseline, which is a rather 
optimistic assumption at first. Fig. [3] demonstrates that 
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Combination with L 


= 3 000 km: 


0.1 


-0.5/+0.5 -1.4/+1.4 


-1.0/+0.9 


-3.0/+2.6 


0.01 


-1.8/+1.7 -5.5/+5.0 


-2.4/+2.3 


-7.3/+6.9 


0.001 


-8.3/+6.9 -27/+21 


-7.9/+7.7 


-26/+25 




Core crossing baseline alone: 




0.1 


-0.5/+0.5 -1.4/+1.4 


-1.1/+1.1 


-3.7/+3.3 


0.01 


-2.1/+5.8 -7.2/+9.2 


-2.5/+2.6 


-7.5/+S.2 


0.001 


-9.9/+19 -40/+35 


-9.5/+11 


-31/+47 



TABLE I: Per cent errors (la and 3a) on pic for an inner core 
crossing baseline and pc for an inner core touching baseline 
as labeled in the columns. The upper group of numbers refers 
to the combination with a baseline L — 3 000 km, the lower 
group to the core crossing baseline alone. The numbers are 
given for different values of sin 2 26n and Sep = 0. 



sin 2 2#i3 > 10~ 3 is a prerequisite to obtain high preci- 
sions, which also applies to the measurement of pc (not 
shown). Since the matter density information mainly 
comes from the appearance channels (Pep and Pep,), it is 
not very surprising that the precision of the measurement 
is suppressed by sin 2 26*i 3 . By comparison of the shaded 
contours (includes correlations) and the dashed curves 
(no correlations), it also demonstrates that the combi- 
nation with L = 3 000 km is very close to the optimal 
(correlation-free) performance. In Tab. [i] the per cent 
errors are listed for different selected values of sin 2 26*i3, 
where also the values for the core crossing baseline alone 
and for the measurement of pc for an inner core touch- 
ing baseline are shown. From Fig. 01 and Tab. [I] one can 
easily see that each core crossing baseline alone is almost 
correlation-free for large values of sin 2 26*i3, whereas for 
small values of sin 2 2#i3 the L = 3 000 km baseline pro- 
vides valuable additional information to reduce correla- 
tions. The reason for this is that the correlation with Sqp 
is rather unimportant for large values of sin 2 2#i3 at these 
very long baselines, where the contribution from the so- 
lar oscillations becomes MSW effect suppressed (see, e.g., 
Ref. 14]). From Tab.|l| we can finally read off a la error 
of less than one per cent for pic under optimal condi- 
tions, i.e., for large sin 2 2#i3 and long enough baselines. 

Let us now restrict these optimal assumptions some- 
what. As it is obvious from Fig.Q] there may not always 
be potential detector locations especially for inner core 
crossing baselines. Therefore, we show in Fig. 0] the pre- 
cision of the measurement for a somewhat smaller value 
of sin 2 26*13 = 0.01 as function of the actual detector lo- 
cation for the laboratories as given in the plot labels. In 
this case, the optimal precision for the inner core crossing 
baselines of about 2% can not be reached because there 
is nothing besides water on the opposite sides of the main 
laboratories. However, a precision of 3% is easily reach- 
able, since, because of the spherical geometry, baselines 
through the inner core travel long distances within the 
inner core for a large region of the inner core shadow re- 
gion projected onto the Earth's surface. The same, in 
principle, applies to the outer core crossing baselines. 

In summary, a precision determination of the Earth's 
absolute inner or outer core density at the percent level 
seems to be feasible with a neutrino factory baseline for 
large enough sin 2 2#i 3 > 10 -2 , even if one takes into ac- 
count the knowledge on the neutrino oscillation parame- 
ters. Since a neutrino factory muon decay ring naturally 
spans two baselines, such a core crossing baseline could 
be an extra payoff of a neutrino factory in addition to its 
main purpose to measure the neutrino oscillation param- 
eters precisely. Alternatively, the optional subsequent op- 
eration of a very long baseline, if more efficient, requires 
the flexibility to modify the storage ring after its main 
purpose has been carried out. Since the major challenge 
is a decay ring with a vertical decay tunnel, a detailed 
feasibility study of these options is needed. From the 
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FIG. 4: (Color online) Precision of the measurement of pc (light gray shaded outer core crossing regions) or pf c (dark 
gray shaded inner core crossing regions) as function of the actual position of the detector for the three laboratories in Fig. 
as neutrino sources (from west to east). The la errors are given within the rings as they are labeled. For all setups, the 
combination with L = 3 000 km, as well as sin 2 26is = 0.01 are assumed. 



geographical point of view, given some of the current po- 
tential candidates for a neutrino factory, there are many 
potential detector locations on the southern hemisphere 
where sufficient precisions could be obtained. From the 
neutrino physics point of view, other applications of such 
a very long baseline would be the degeneracy-free mea- 
surement of sin 2 26*13 j the direct verification of the MSW 
effect, or the parametric enhancement of neutrino oscil- 
lations [lij . 

In comparison to seismic wave geophysics, neutrinos 
are sensitive to the absolute electron density in mat- 
ter, which means that the relationship to the absolute 
matter density is much cleaner from model-dependent 
assumptions. This measurement could therefore help 
to test the equation of state for seismic waves. The 
obtainable precision at the percent level is competitive 
to the precisions of the density jumps given by seismic 
wave geophysics. For example, at the inner core bound- 
ary, Ap ~ 0.82 ± 0.18 gem -3 0. Of course, the neu- 
trino factory approach alone is unlikely to give a three- 
dimensional model of the Earth, but its strength to mea- 
sure the baseline averaged density as opposed to the vol- 
ume averaged density makes it a good candidate for di- 
rect matter density tests of especially the innermost parts 
of the Earth. 
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